Laser control of complete vibrational transfer in Na 2 using resonance coalescence 
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With a specific choice of laser parameters resulting into a so-called exceptional point in the 
wavelength-intensity plane, it is possible to produce the coalescence of two Floquet resonances 
describing the photodissociation of the molecule Na2, which is one of the candidates for molecular 
cooling. Appropriately tuning laser parameters, following a contour around the exceptional point, 
the resonances exchange their labels. This represents a laser control of the vibrational transfer 
from one field-free state to another, through an adiabatic transport involving these resonances. The 
proportion of undissociated molecules at the end of the pulse is checked through Floquet adiabatic 
theory. A vibrational cooling scenario can be proposed based on a complete vibrational transfer 
which is predicted, with only 20 percent of molecules undergoing dissociation. 

PACS numbers: 33.80.-b, 42.50.Hz, 33.80.Gj 
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Progress made during the last two decades in the cool- 
ing of quantum gases at ultracold temperatures down to 
the quantum degenerate regime has opened a possibil- 
ity to study small molecules at the level of individual 
quantum rovibrational and hyperfine states [3, 0|. An 
important step of the experimental developments is to 
achieve the deterministic manipulation of quantum states 
of small molecules, the formation of pre-selected quan- 
tum states and the control of the coupling between such 
states. The variety of robust experimental techniques 
[H, |i[ which are currently being developed may later be 
used as a tool set for many-body physics, high preci- 
sion measurements and technological applications such 
as quantum computing. Ultracold alkali-metal molecules 
are most often prepared in a distribution of vibrational 
levels of the singlet electronic ground state, or of the 
lowest triplet state. An important challenge is to cre- 
ate dense samples of ultracold molecules in a selected 
single quantum state. At present, there are many dif- 
ferent techniques used to create and manipulate well- 
defined rovibrational and/or hyperfine states of diatomic 
molecules at ultralow temperatures, e.g. photoassocia- 
tion (PA) Jp, stimulated Raman adiabatic passage (STI- 
RAP) [6HaK black-body radiation-assisted laser cooling 
(TSELor optical pumping by a broadband femtosecond 
laserpj,[l3|. 

In this work we propose an alternative way for a de- 
terministic manipulation of vibrational states of alkali- 
metal diatomic molecules, relying on adiabatic transfer 
of population using chirped picosecond laser pulses, blue- 
detuned with respect to an atomic transition frequency. 
The technique is based on the coalescence of two res- 
onance energies described using the Floquet formalism 
1141. The laser-controlled vibrational transfer scenario 
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consists in applying an electromagnetic field which cou- 
ples the ground or a metastable electronic molecular state 
to an excited electronic state. The induced resonances 
are identified through a series of complex eigenenergies of 
the Hamiltonian of the system, the imaginary part being 
related to their dissociation lifetime. With an appropri- 
ate choice of the laser parameters (intensity, wavelength) 
which is called an exceptional point (EP) [l^, [l6[ of the 
parameter space, it is possible to produce a degen- 
eracy of two such energies. A specific loop encircling the 
EP in the laser parameter plane has to be identified, such 
that starting from a given field-free vibrational state, a 
single resonance adjusts continuously its characteristics 
(energy and width) to reach a different field-free vibra- 
tional state [lH • In a second step we introduce the 
time requested to follow the loop and fulfill the condi- 
tions for an adiabatic transfer. This finally leads to an 
estimate of the percentage of undissociated molecules at 
the end of the pulse. 

The proposed scheme is general but we consider here- 
after the example of the Na2 molecule in an excited vi- 
brational level of its a 3 S+ lowest triplet electronic state 
(referred to as state u), which can be formed in optical 
traps |20J . The formation of ultracold sodium molecules 
in particular, has been demonstrated in several experi- 
ments |2l| - |23| , and therefore the present application rep- 
resents a realistic situation. The same technique can be 
applied to other alkali dimers formed in optical traps with 
only small changes due to the similarities of molecular 
potentials for all alkali dimers. The laser-controlled vi- 
brational transfer scenario consists in applying an elec- 
tromagnetic field with wavelength around 560 nm that 
couples state u with the (l) 3 H g (3 2 S + 3 2 P) excited elec- 
tronic state, labeled g below. As a complement to our 
preliminary study on (X 2 £+ — > A 2 E+ photodissoci- 
ation), the present situation looks quite promising. The 
reduced mass is about 20 times larger and the poten- 



tial well of the initial state is about 200 times shallower, 
than in the case, so that the density of vibrational 
levels is about 2 orders of magnitude larger than in H^. 
We thus expect resonance coalescence to occur at signif- 
icantly smaller intensities, while the proximity of levels 
will require larger pulse times to remain in the adiabatic 
regime. Ultimately, we show that about 80% of the Na2 
molecules survive to photodissociation, in contrast with 
the 10% amount found for Hj. 

Hereafter we refer to a model describing a rotationless 
field-aligned molecule in one spatial dimension (the in- 
ternuclear distance R) and involving only two electronic 
states \u) and \g) whose Born-Oppenheimer (BO) poten- 
tial energy curves V® ° (R) and V g BO (R) are displayed in 
Fig. QJ,, b. Among all electronic states of Na2, the cho- 
sen \u) and \g) states present a favorable curve crossing 
situation at the laser wavelengths used to couple them 
(Fig. [TJ:). In the following we neglect the spin-orbit in- 
teraction in the Na2 molecule as it does not change the 
overall picture. The frozen rotation assumption can be 
validated by considering the long rotational periods of 
Na2 (estimated as hundreds of ps) as compared with the 
pulse durations under consideration (less than 2 ps) . The 
time-dependent wave function of the system is expanded 
on these two states 



\y(R,t))= X u(R,t)\u)+ Xg (Rit)\g) 



(1) 



with the unknown functions \u, Xg accounting for the 
field- assisted nuclear dynamics. When the interaction 
of the molecule with a continuous wave (cw) field with 
pulsation lo is described by a time-periodic Hamiltonian, 
applying the Floquet ansatz, the two-component solution 
of the time-dependent Schrodinger equation (TDSE) can 
be written as 



(2) 



where Ep is the Floquet quasi-energy. The periodicity 
in time of 4>k{R,t) (k — u,g) allows us to expand these 
functions in Fourier series 
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Finally, adopting the length-gauge and the long- 
wavelength approximation for the matter-field coupling, 
the Fourier components ip 7 ^(R) are given as solutions of 
the set of coupled differential equations resulting from 
the TDSE, exemplified by the general form, for all n [l4| 
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where T = (~h 2 /(2M))(d 2 /dR 2 ) is the usual radial ki- 
netic energy operator and A4 is the reduced mass. The 
transition between the u and g states is governed by the 
independent molecular transition dipole moment /«(-R) 
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FIG. 1. (Color online) The Na2 potential curves involved in 
the proposed scheme, (a) the a 3 E^(3 2 5' + 3 2 S), and (b) the 
(1) 3 H g (3 2 S + 3 2 P) potentials. The absorption of a photon 
with energy hw — 0.0808 a.u. is pictured with a vertical 
arrow, (c): the same potentials dressed by the laser field to 
produce (g, -1) and (u, 0) channels. The adiabatic potentials 
(black lines) obtained by a diagonalization of the potential 
matrix including the laser-induced coupling are also shown. 
Horizontal lines indicate positions of bound and quasi-bound 
levels obtained after solving the Schrodinger equation for the 
corresponding one-channel adiabatic (black) or diabatic (red) 
potential curves. 



[2J, [25j. The laser is characterized by the linearly- 
polarized electric field £(t) = £ n cos(ujt), the wave- 
length A = 2itc/uj and the intensity I oc £q. Solu- 
tions with Siegert outgoing wave boundary conditions 
in the open channels produce complex quasienergies of 
the form Ep = Er — iTn/2 where Tr is the resonance 
width related to its decay rate. The single photon pro- 
cesses are described using the two diabatic channels u 
and g, with field-dressed potentials V U (R) = V^°(R), 
V g (R) ee V g BO {R) - hw (FigQJ:), and with wave func- 
tions fu(R) an d fg 1 (-R), designated briefly as <p u (R) and 
f g {R)- It is of course possible, for each value of the coor- 
dinate R, to obtain, through a linear transformation, the 
functions associated with the adiabatic potentials V- (R) , 
V+(R) (FigOJ) resulting from the diagonalization of the 
diabatic potential matrix after accounting for the matter- 
field interaction (wavelength A = 563.903 nm and inten- 
sity / = 0.35 GW/cm 2 ) H). 

The adiabatic potentials are very useful for the inter- 
pretation of resonances. In particular, for high laser in- 
tensities we expect the states associated with these adi- 
abatic potentials to be very good zeroth-order approxi- 
mations for the resonances of the full problem. We have 
previously observed [27j that the Floquet resonances in- 
duced by high intensity electromagnetic fields fall into 
two categories: (i) the Feshbach-type (FT) resonances as- 
sociated with the levels of the V+ (R) potential well inter- 



acting with the continuum of the V_ (R) potential; (ii) the 
shape-type (ST) resonances, decaying through or above 
the barrier of V- (R) and still interacting with V+ (R) by 
residual non-adiabatic kinetic couplings. The energy of 
FT resonances increases with increasing intensity, as the 
adiabatic potentials split further and further; these res- 
onances ultimately merge with the vibrational levels of 
V+ (R) of vanishing width for very high intensity. Simi- 
larly the energy of ST resonances decreases with increas- 
ing intensity as the barrier height of the V_ (R) potential 
decreases. Their widths grow until the flattened barrier 
is no more able to accommodate a resonance. 

In order to locate the EP resulting from the coales- 
cence of FT and ST resonances, we calculate the reso- 
nance energies as a function of intensity for various val- 
ues of the wavelength [13, [3 ■ Our best estimate for the 
parameters leading to coalescence of the v = 3, 4 pair is 
X E p = 562.53 nm and I E p = 0.332 GW/cm 2 . In con- 
formity with our previous analysis (l7j . for A < Xep the 
widths as a function of intensity, exhibit an avoided cross- 
ing with a tweezer-like profile, while the energies cross 
each other at the EP. These patterns are exchanged be- 
tween widths and energies for A > Xep ■ This means also 
that with a small change of the wavelength (a slightly 
larger value in our case) we can produce an avoided 
crossing for the energies while the rates will be cross- 
ing. According to it should be possible to find a 
contour in the laser parameter plane that goes continu- 
ously from one resonance to another around an EP. A 
workable control scheme would be to apply adiabatically 
a chirped pulse with the goal to transfer population from 
one field-free state to another. Adiabaticity refers here 
to a strategy consisting in the transport of a vibrational 
level on a single Floquet resonance, following it in time 
until it changes its label and ends up into another single 
vibrational level. This requires a laser pulse duration tf 
that exceeds the vibrational transition time, roughly es- 
timated as h/AE v , AE V being the energy gap between 
adjacent vibrational levels. The resulting penalty, when 
taking a long duration laser pulse, is a possible deple- 
tion of molecular bound states. It is thus very important 
to reach a compromise between adiabaticity and the re- 
quirement for a dissociation not to be complete along the 
contour. We have built such a contour with an intensity 
and a wavelength obeying to a simple relation 



I = I» 



sin(</>/2), A = A + 5\ sin(^) , (5) 



with Ao = 562.53 nm, i.e. the wavelength for the EP 
of the pair 3,4, I max = 0.4 GW/cm 2 and 8X = 2 nm. 
The parameter <fi is continuously and linearly changing 
with time from (at t = 0) to 2tt (at t = tf). It is 
important to emphasize the experimental feasibility of 
such a picosecond laser pulse in the visible domain, with 
a chirped amplitude around 560 nm not exceeding 1% 
and a rather modest intensity, less than 0.5 GW/cm 2 . 
Figure [5] shows the loop in the parameter plane (panel a) 
and the resulting trajectories in the energy plane (panel 
b). The trajectories going from v = 3 to v — 4 or from 
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FIG. 2. (Color online) (a) The laser pulse in the parame- 
ter plane (A, /) producing (b) the trajectories of the energy 
plane either starting from the field-free vibrational level v = 3 
(black dashed line) or v = 4 (red solid line). Because the loop 
goes around the EP, there is a resonance exchange, (c) The 
probability Pnd for the molecule to survive after being ex- 
posed to 2 chirped pulses of duration 835 (dashed lines) and 
1670 femtoseconds (solid lines). Two cases are considered: 
transitions from v = 3 to v = 4 (black lines) and from v = 4 
to v — 3 (red lines). 



v = 4 to v — 3, involve different rates, in agreement with 
the previous analysis of ST and FT resonances 17 1 . 

In order to estimate the fraction of undissociated 
molecules left after the pulse is over, we refer to the adi- 
abatic Floquet theory [27], [Hj]. This assumes that the 



chirped laser pulse envelope and its frequency vary suffi- 
ciently slowly with time such that the overall fraction of 
nondissociated molecules -Pnd is given by 



Pnd(*) = exp 



T R (t') dt' 



(6) 



Here, r#(t) is associated with the relevant Floquet quasi- 
energy eigenvalue calculated using the instantaneous field 
parameters at time t. We present in Fig. [2J: Pnd for two 
cases: either starting from v — 3 to reach v — 4, or the 
reverse process and for two different laser durations tf. 
The results differ drastically. As already mentioned, the 
interpretation relies on the nature of the instantaneous 
resonance states created during the pulse. In the first 
case (v = 3 to v = 4), during the first half of the pulse 
with increasing intensity and A > Xep the resonance is 
ST, with large values of the decay rates. In the second 
half of the pulse with decreasing intensity it also matches 
with the ST resonance v = 3. For the transfer v = 4 to 
v = 3 the situation is opposite: The FT resonance of the 
first half of the pulse matches with the FT resonance of 
the second half, after the exchange of label when going 
across Xep- Therefore there is a better chance to pre- 
serve the molecules from dissociating when the transfer 
amounts to a decrease of the vibrational quantum num- 
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ber. 

A crucial result of this work is that at the end of 
the pulse more than 80% of the molecules survive the 
dissociation and are completely transferred from vibra- 
tional state v = 4 to v = 3. The importance of this 
result can be understood by comparing with a similar 
vibrational transfer control in the lighter molecular sys- 
tem H^. As previously mentioned, due to much higher 
density of levels, resonance coalescences in Na2 are ob- 
tained for much lower laser intensities. More specifically, 
I EP ~ 0.3 GW/cm 2 for Na 2 , as compared with I EP ~ 0.4 
TW/cm 2 for H^. As a result the resonances participat- 
ing in the adiabatic transfer process have decay rates of 
the order of ten cm -1 while, in the case of H^~, they 
reach a thousand of cm -1 , leading thus to much better 
protection against dissociation for Na 2 . However, a high 
density of levels requires longer pulse durations, for the 
adiabaticity to be safely implemented. More precisely, 
50 fs pulse was long enough for an adiabatic transfer in 
, whereas about 1000 fs is necessary for the same adi- 
abaticity to be achieved in Na 2 , which is thus submitted 
to laser interaction and dissociation for a duration about 
20 times longer than H^. Despite this fact, the long 
lifetimes of Na 2 Feshbach resonances on which the vibra- 
tional transfer relies are such that 80% of Na 2 molecules 
remain non-dissociated at the end of the pulse, as com- 
pared to only 10% in the case of H^. It is worth noticing 



that there are other lifetime- limiting processes competing 
with the photodissociation, namely spontaneous emission 
(for the (1) 3 II 9 state) and the dimer-atom inelastic colli- 
sions between Na and Na 2 . The Na 2 lifetimes associated 
with these processes are much larger than 1000 fs, of the 
order of a few ns for the spontaneous emission and a few 
seconds for inelastic atom-dimer collisions. 

We have shown that complete vibrational transfer is 
possible using resonance coalescence (exceptional points) 
in laser-induced photodissociation/photoassociation of a 
molecular system using a rather simple experimental 
technique. This is illustrated for the case of Na 2 with re- 
alistic laser parameters. Comparing the cases of and 
Na 2 , we have shown that the proposed scheme can be ap- 
plied to other molecular species like Cs 2 or Rb 2 , because 
the wavelength and the amplitude of the applied electric 
field are two parameters giving considerable freedom to 
manipulate the resonance quasienergies. The ultimate 
goal for this and future similar studies is molecular vi- 
brational cooling using a purification scheme to obtain a 
single ro-vibrational state, by successive transfers start- 
ing, for instance, from a thermal distribution of initial 
vibrational levels. 
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